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Abstract 

We present an algorithm to solve BSDEs with jumps based on Wiener Chaos Expansion and 
Picard’s iterations. This paper extends the results given in [B] to the case of BSDEs with 
jumps. We get a forward scheme where the conditional expectations are easily computed 
thanks to chaos decomposition formulas. Concerning the error, we derive explicit bounds 
with respect to the number of chaos, the discretization time step and the number of Monte 
Carlo simulations. We also present numerical experiments. We obtain very encouraging 
results in terms of speed and accuracy. 
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1. Introduction 

In this paper we are interested in the numerical approximation of solutions {Y, Z, U) to 
backward stochastic differential equations (BSDEs in the sequel) with jumps of the following 
form 

Yt = i+ r f{s, w, Us) ds - r ZsdBs - 

Jt Jt 

where i? is a 1-dimensional standard Brownian motion and A is a compensated Poisson 
process independent from B, i.e. Nt := Nt — nt and {Nt)t is a Poisson process with intensity 
K > 0. The terminal condition is a real-valued J^T^measurable random variable where 
{Bt}o<t<T stands for the augmented natural hltration associated with B and N. Under 
standard Lipschitz assumptions on the driver /, the existence and uniqueness of the solution 
have been stated by Tang and Li |2I], generalizing the seminal paper of Pardoux and Peng 

pE]. 

The main objective of this paper is to propose a numerical method to approximate the 
solution {Y,Z,U) of ([^. In the no-jump case, there exist several methods to simulate 
(U, Z). The most popular one is the method based on the dynamic programming equation, 
introduced by Briand, Delyon and Memin [5]. In the Markovian case, the rate of convergence 
of the method has been studied by Zhang [22] and Bouchard and Touzi [3] . From a numerical 
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point of view, the main difficulty in solving BSDEs is to compute conditional expectations. 
Different approaches have been proposed: Malliavin calculus [3], regression methods P] and 
quantization techniques [I]. In the general case (i.e. for a terminal condition which is 
not necessarily Markovian), Briand and Labart [B] have proposed a forward scheme based 
on Wiener chaos expansion and Picard’s iterations. Thanks to the chaos decomposition 
formulas, conditional expectations are easily computed, which leads to an efficient, fully 
implementable scheme. In case of BSDEs driven by a Poisson random measure, Bouchard 
and Elie P] have proposed a scheme based on the dynamic programming equation and 
studied the rate of convergence of the method when the terminal condition is given by 
^ = g{XT), where 5 ^ is a Lipschitz function and X is a forward process. More recently, 
Geiss and Steinicke |H] have extended this result to the case of a terminal condition which 
may be a Borel function of finitely many increments of the Levy forward process X which 
is not necessarily Lipschitz but only satishes a fractional smoothness condition. In the case 
of jumps driven by a compensated Poisson process, Lejay, Mordecki and Torres 1131 have 
developed a fully implementable scheme based on a random binomial tree, following the 
approach proposed by Briand, Delyon and Memin |1]. 

In this paper, we extend the algorithm based on Picard’s iterations and Wiener chaos 
expansion introduced in [B] to the case of BSDEs with jumps. Our starting point is the use 
of Picard’s iterations: (F°, U^) = (0, 0, 0) and for g e N, 

Y,-^' = (+ r f{B, 1 ?. Zl,U^)ds - r-dB,- [ UI+'dN„ 0 < i < r. 

Jt Jt 

Writing this Picard scheme in a forward way gives 


= EU + J^f ( 5 , y/, zi U^) ds I Xt) - / (s, y/, zi u^) ds, 






Xd- 





Xd: 


+ [f(s,Y;,ziut)ds 
jY(s,Y;,ziui)ds 



where Df'^X (resp. D^^^X) stands for the Malliavin derivative of the random variable X 
with respect to the Brownian motion (resp. w.r.t. the Poisson process). 

In order to compute the previous conditional expectation, we use a Wiener chaos expan¬ 
sion of the random variable 

F^ = i+ r f{s,Y^,ZlU^^)ds. 

Jo 


More precisely, we use the following orthogonal decomposition of the random variable F'^ 
(see Proposition 2.6) 


00 k 


k=i /=o k,eN' 
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where (resp. Ll^"’^{g)) denotes the iterated integral of order m oi g w.r.t. the 

Brownian motion (resp. w.r.t. the compensated Poisson process), (e[A:i,..., /cm])fc,„eN is an 
orthogonal basis of (L^)®'”([0, T]), the subspace of symmetric functions from (L^)®™’([0, T]). 
The sequence of coefficients (dk;,jj._Jk;eN',jfe_;eN'=-' ensues from the Wiener chaos decomposi¬ 
tion of F^. 

The point to get an implementable scheme is that we only keep a finite number of terms 
in this expansion: we use a hnite number of chaos and we choose a hnite number of functions 
(ei, • • ■ , Cat) to build (e[fci, • • • , , 7 V}- More precisely, if we choose e* := 

where = ih and h := we obtain 


N 


F*? ~ E [F^ 


+ EE 

k=\ \n\=k 


di n K. 


B-t -Bi 


2=1 




C^AN.-N.^Kh) 


where Ki (resp. Ci) denotes the Hermite (resp. Charlier) polynomial of degree i, n = 

,n^) is a vector of integers and \n\ = -1- nf). By using this 


B 
1 ) 




approximation of we can easily compute E(F'^|Ft), E(F{°^F'^|Ft-) and E(F)^''F'?|Ft-), 
which gives us ^ . To get a fully implementable algorithm, it remains to 

approximate E(F'^) and the coefficients {dAn,k by Monte Carlo. 

When extending |B] to the jump case one realizes that the main difficulty lies in the fact 
that there is no hypercontractivity property in the Poisson chaos decomposition case. This 
property plays an important role in the proof of the convergence in the Brownian case. To 
circumvent this problem, we exploit a recent result of Last, Penrose, Schulte and Thale nn, 
which gives a formula to compute the expectation of products of Poisson multiple integrals, 
and the according result for the Brownian case from Peccati and Taqqu im. In fact, in 


i(i) 


equation (16) of Proposition 2.9 we get an explicit expression for 


(/„,)) 


in terms of a combinatoric sum of tensor products of the chaos kernels /„.. Here Imifni) 
denotes the multiple integral of order rii with respect to the process B + N. By this expres¬ 
sion one gets the required estimates for the truncated chaos without the hypercontractivity 
property. Therefore, to prove the convergence of the method we may proceed similarly to 
[0], and split the error into four terms: 

• the error due to Picard iterations 


• the error due to the truncation onto the chaos up to order p 

• the error due to the hnite number of basis functions (ci, • • ■ , ca?) for each chaos 

• the error due to the Monte Carlo simulations to approximate the expectations appear¬ 
ing in the coefficients {d1)n,k- 


The paper is organized as follows: Section contains the notations and gives preliminary 
results. Section describes the approximation procedure. Section states the convergence 
results and Section presents the algorithm and some numerical examples. Some technical 
results are proved in the appendix. 
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1.1. Definitions and Notations 

Given a probability space (G, P) we consider 

• := P); P £ bl*; th® space of all J^r-measurable random variables (r.v. 

in the following) X : G i—)■ M satisfying ||X||p := E(|X|^) < oo. 

• S^(M), p E N,p > 2, the space of all cadlag, adapted processes 0 : G x [0,T] i —> M 
such that || 0 ||gp = E(sup^g[o,r] \4>t\^) < C) 0 - 

• p E N,p > 2, the space of all predictable processes 0 : G x [0,T] i —> M such 
that || 0 ||hp = E(/(J’ | 0 i|P(if) < oo. 

• L^(0,T), the space of all square integrable functions on [0,T]. 

• C^’\ the set of continuously differentiable functions 0 : {t,x) E [0,T] x with con¬ 
tinuous derivatives w.r.t. t (resp. w.r.t. x) up to order k (resp. up to order 1). 

• C^’\ the set of continuously differentiable functions 0 : (t, x) E [0, T] x with contin¬ 
uous and uniformly bounded derivatives w.r.t. t (resp. w.r.t. x) up to order k (resp. 
up to order 1). The function 0 is also bounded. 

• Il^sp/IIL) the sum of the squared norms of the derivatives of /([0,T] x M^,M) w.r.t. all 

the space variables x which sum equals j : || 0 fp/||^ := J2\k\=j where 

\k\ = ki + k2 + fcs- 

• Cff, the set of smooth functions / : M"' i —> M (n > 1) with partial derivatives of 
polynomial growth. 

• IK') ■) OIIlp) P ^ norm on the space Sy(M) x H^(M) x H^(M) dehned by 

\\{Y,Z,U)\\l, :=E( sup \Ytn+ E{\Zt\ndt + k E{\Ut\ndt. (2) 

te[o,T] “'0 Jo 


Hypothesis 1.1. ITe assume 

• the terminal condition ^ belongs to 

• the generator f E G([0,T] x M^;M) is Lipschitz continuous in space, uniformly in t: 
there exists a constant Lf such that 

\f{t,yi,zi,ui) - f{t,y2,Z2,U2)\ < Lf {\yi - yal + - ^2! + \ui - M2I) • 


Lemma 1.2. If Hypothesis 1.1 is satisfied and ^ E (defined below) we get from 
Theorem S.f] that for a.e. t E [0,T] 


Zt = E[D^PYt\Ht-l Ut = E[DY%\Ht-] P-a.s 


( 1 ). 


(3) 


where Df'^X stands for the Malliavin derivative w.r.t. the Brownian motion of the random 
variable X, and D^^^X stands for the Malliavin derivative w.r.t. the Poisson process of the 
random variable X. Here E[-|J^t_] should be understood as the predictable projection, and 
since the paths s 1 —>■ D^'^Yg are a.s. cadlag we define D^f^Yt := hms 4 .t if the limit exists, 

and zero otherwise. 
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2. Wiener Chaos Expansion 

2.1. Notations and useful results 

2.1.1. Iterated integrals 

We refer to m and 1121 for more details on this section. Let us briefly recall the Wiener 
chaos expansion in the case of a real-valued Brownian motion and an independent Poisson 
process with intensity k > 0 . 

We dehne 


Go(t) — Bt, Gi(t) — Nt — K,t, 


and L*!’"'’*'“(/) the iterated integral of / with respect to Go and Gi 

= If [If ■■■[If /(ii,....4MG„(0)V"dG4_,(4-i)) dGi^ih). 

We have the following chaotic representation property 
Proposition 2.1. (fm Proposition 2.1]) Fork E N* define 

ifc := {ii,...,ik) e {0,1}^ 

Any F E L‘^{Ft) has a unique representation of the form 

OO 

F=E(F) + i: Lh/ij (4) 

k=i ifeeto,!}* 

where E and = {(ti, ... fik) G [0,T]^ : 0 < < ■ ■ ■ < tfc < T} is the simplex 

o/[0,T]^ 

Let I ifc I := Due to the isometry property it holds 

llL‘(/)f = '="‘'11/111.. 

and for any / G Lfi{fLk)i 9 ^ ifc ^ { 0 ; 1 }^) and jm G { 0 , 1 }™ we have (see [Ml 

Proposition 1.1]) 


ElL'{{f)Lt(g)] = { 

Then, ||Fp = E[F]" + Ej>i 
p is dehned by 


/sfe fiG, • • • , tk)g{ti, ■ ■ • , tk)dti ■■■dtk if ifc = jm 

0 otherwise. 

Dhe chaos approximation of F up to order 


4(F) :=E(F) + 1:J:4>(/,.). 

k=l i/g 


(5) 


We also dehne Pfc(F) := Ei, Ll^ifh)- We have 

E[('Pfc(F))^]=E'^'‘'''ll/iJI| 

ifc 


( 6 ) 
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• Let / G L^{T,k) and j G {0,1}. Following [T3], we define the derivative of w.r.t. 

the Brownian motion and the Poisson process as the element of L^{Q x [0,T]) given 
by 




1=1 


1-1 


where i means that the i-th index is omitted. 

• Let j G {0,1}. We extend the dehnition of to 

{ oo k 

F G L‘^{Ft) satisfying @ and ll/iJIlfe < 

k=l ife /=1 

If F G Dom then 

- E|Fp + -t'E r\Dl»F\^dt < oo. 

^ 0 

• F with chaotic representation Q belongs to Dom D =: if F belongs to 

Dom D^^^ n Dom D^^\ i.e. 

OO 

||F|li... := E|Fp + E *:E''"‘'ll/iJII. < 

k=l i/c 

More generally, we dehne as follows: 

• Let m > 1. We say that F satisfying Q belongs to if it holds 

m oo 7 I 

||F||^„,, := E|Fp + EE Tr^ni E ''"‘'ll AJIL < «>■ 

1=1 k=l ^ ^ * ifc 

We recall 

We dehne for I G N* the seminorm II ■ \\jji on by 


IIFiri,. :=E''"‘'E 




■ 


dti - ■■dti] =Y^ 


l|2 


(7) 


where Dl[ ... j-^ = Dl^ ■ ■ ■ Dl\ represents the mnlti-index Malliavin derivative. 
Remark 2.2. By using this notation we have ||F ||^„,,2 = ]E|Fp + YaLi ||-^|li)i- 
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• For m > 1 and j G N* we define D™-’-? as the space of all F G D"*’^ snch that 

\\^\\L,j-= esssup < cx) 

l<l<m iie{0,l}* {tl,- ,tl)&[0,TY 


^rn,j (denotes the space of all triples of processes {Y, Z, U) belonging to S:p(M) x (M'^) x 


H^(]R) and snch that 


\\{Y,z,u)\\i„y.= •£ •£ esssup < «. 

i<i<m i; hi,-” ,ti)e[o,r]* 


where || ■ \\\^j has been dehned in ([^. We denote 5™’°° ;= 

Remark 2.3. If F := g{G), where : M —)■ M zs a function and G G D^’^, we have 
(following Proposition 5.1]) that 

{Di^^F, Di^^F) = {g'{G)Df^G, g{G + D^^^G) - g{G)). 

Moreover, using Notation Q, we get 

IIFIIP = + K||£>">F|li,,u>,|„,Ti) < IIS'IILIIGIIL.. 

More generally, if g : ^ M. is a CJf' function and G G we have 


D” 


Lemma 2.4. Let 1 < m < p + 1 and F G D™’^. We have 

E[\F - C,{F)\^] < 

Proof. Using ([^, we get 


w^r 


(p + 2 - m) ■ ■ ■ (p + 1)' 


E[|f-A(An= i: Ein(F)"]= i: 


ifc 11 


k>p-\-l fc>p+l ife 

k\ [k — m)\ 


„^_,,{k-m)\ k\ 


< 


k\ 


< 


(p + 2 - m) ■ ■ ■ (p + 1) 
1 


E 




k>p+l V Z ij, 

hi 

V V/-rl‘fel|l f- l |2 

(p + 2 -m)---(p+l),t;(h-m)!^ 


A 


□ 
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2.1.2. Multiple integrals 

In the following, A denotes the Lebesgue measure. Setting 

M{ds,dx) := dGo{s)d6o{x) + dGi{s)dSi{x) 

we get an independent random measure in the sense of ltd (see HU]). There exists a chaotic 
representation by multiple integrals w.r.t. this random measure M which is equivalent to 


Proposition 2.1 


Proposition 2.5. ffTOf ) Any F G L^{Ft) can be represented as 

CO 

F = E[F] + Y,higk). 


k=l 


With Qk e (L2)®fc(A ® (5o + K(5i)) := (L2)®^([o,T] x {0,1}, S([0,T]) ® A ® (5o + k^i)). 
This representation is unigue if we assume that the functions gk{zi, ...,Zk) with Zi = {ti,Xi) G 
10, n X {0,1} are symmetric. 


In fact it holds for symmetric 

hidk) = - ,(-,4))), 


( 8 ) 


ife 


where is dehned in Proposition 2.1, and for symmetric g^ G (L^)®*^(A ® ((5o + k4)) and 
U e (L2)®-(a 0 (5o + «5i)) 




If F G then we have 


gk{{ti, ii),..., {tk, ik)) = k<m. 


(9) 


( 10 ) 


For the implementation of the numerical scheme we will use Hermite and Charlier poly¬ 
nomials. In order to do so, we provide a chaotic representation consisting only of iterated 


integrals of the form and for which the relations (20) and (21) below can be 

used. 

Use {po,pi} = {l{o}; as orthonormal basis of T^({0,1}, 2^°’^^, ho + ^^i) and £x an 

orthonormal basis {ek)keN for L^([0, T], i3([0, T]), A). By setting 

e[{ki,ii), ■ ■ ■, {km,im)] ■= (e^i ^Pn) ® ® (e^^ <S)Pi^), kj G N,ij e {0,1} 

we get an orthonormal basis of (L^)®™'(A ® (ho + ^hi)). The symmetrizations 

e[(/Ci, il), . . . ) (^m; ^m)] • ^ | ^ ) ®[(^7r(l)) 4(1)); ■ ■ ■ i (^7r(m)) 4(m))]) hj G N, ij G {0, 1} 

( 11 ) 


TT&Sn 






0 ?^ 


form an orthogonal basis of (L^) (A 0 (^o + the subspace of symmetric functions from 
(L2)®-(A0((5o + K(5i)). 

We also will use the notation 


l[ki, . . . , km] ■— , ^ 


ml 


TV^Sm 


( 1 ) 


Cfc , ,, kj e N. 


Proposition 2.6. Any F G L‘^{Ft) can be represented as 


oo k 


F^EIF] + J2Y1 H H (e|ji. 

k=ii=o ki&W jfe.jeNfc-' 


Jk-l])- 


where 




Proof. According to cm Theorem 1] a permutation of the coordinates of the kernels does 
not change the multiple integral, i.e. for any tt G we have 


^kic\{k\.i A), • • • , (fcfc, ^fe)]) (c[(^7r(l)) ^7r(l))) • • • ) ik'F{k)i ^7r(fc))])- 

For any vr with (v(i),..., iT,(k)) = (0, • • •, 0,1,..., 1) (we assume that (A, ■ ■ ■ ,ik) contains I 


zeros) it holds by the product formula for multiple integrals ( see the Appendix Appendix 


A.5 or cm Theorem 3.6]) 


^l)) ■ ■ ■ 1 {kk^ ^fc)]) dll{e\fk.j^(\p 0 ), . . . , (fc7r(i), 0 )])/fc—;(e[(fc 7 r(/_|_ip 1 ), . . . , (fc7r(fc), 1 )]) 

( 12 ) 

since 

^[(^7r(l): ^7r(l))5 * * * 5 (^7r(fc)5 ^7r(/c))] ^[(^7r(l)7 b): • • • 7 (^7r(/)7 b)] ^ ^[(^7r(/+l)7 I): • • • 7 (.^n{k)^ ^)]7 


and for the contraction-identihcation 0 ])^ (for the dehnition see (A. 12 )) it holds 

€'[(^7r(l)7 b)7 ■ ■ • 7 {k-n{l)-i b)] ®m ^[(^Tri^+l)7 1)7 • • • 7 (^7r(fc)7 1)] b 

if r 7 ^ 0 or m 7 ^ 0. Since 

c[(^ 7 r(i+l )7 1)7 • • • 7 {f^n{k)y 1 )] = k-l c[kT^g^i), . . . , kj^-l^k)] 

K 2 


we conclude from ( 12 ) and ([^ that 

l\(h — ly ^ ^ ^ ^ 

.f^fc(c[(^l7 ^ 1)7 • • • 7 {^ki Ik)]) feVi (c[^7r(l)7 • • • 7 ^■ 7 r(p]).h^_^ (e[/c.n-(;_|_ip . . . , A;7r(fc)]). 


k — l 

K 2 


(13) 


The symmetric functions pk from Proposition |2.5 can be written as 

k 

9k EEE {hki c [{ ki , b), [ ki , 0 )] ® e[(jl 7 1)7 ■••7 { jk-h l)]){L2)®fc 

* = 0 k; jk_i 

X e[(/ci, 0 ), { ki , 0 ), (ji, 1 ), { jk-u l)]ck,,jfc_, 


(14) 
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where we sum over all k; G N* and G ^ and 

— ||c[(fcl, 0), ..., (fc;, 0), (jl, 1), ..., {jk-h 1)] II 

denotes the normalizing factor. 

Lemma 2.7. Fix N eN* and let 

N N 

e[{ki, 0), {k, 0), (ji, 1), {jk-i, 1)] = 0(ei ® ® ® 

i=l j=l 

i.e. nf and nf (1 < ^ < iV) denote the multiplicities of the functions Cj ® Po o.'^d Ci ®pi, 
respectively, so that \n^\ = I and \n^\ = k — l. Letn^l := n^\ ■ ■ ■n^! for A = B,P and define 
n := {n^,n^) so that \n\ = \n^\ + |n'^|. Then 

2 1^1* 

~ Il^[(^l5 0)7 •••5 0)5 1)5 •••5 l)]il(L2)®fc(A(g)(5o+K5i)) ~ ^P] * 

Proof. To compute Ck;,jj._; notice that the functions hj := (cfc^. ® Pif) and hp (1 < j, j' < k) 
are either equal or orthogonal in Lf{\ ® (hg + i^di)). Denoting 

e[{ki, 0 ), {ki, 0 ), (ji, 1 ), {jk-i, 1 )] = hi® ■ ■ ■ ® hk 


yields 


\e[{ki, 0), {ki, 0), (ji, 1), {jk-h l)]II^L2)®'=(A(gi(<5o+K5i)) 


n^\n^\ ^ 


h. 


7v{k) 


(L2)®fc(A®(,5o+K5i)) 


k\ 

n^\n^\ 

k\ 


N 


(8)(d ® Po)®”* ® (8)(ei ® Po)®”^' 

i=l j=l 


(L2)®fe(A®(5o+K<5i)) 


□ 


Abbreviating 


d 




l\{k-l)\ 

K 2 


{9k, e[(/ci, 0 ),{ku 0)] ® e[(ji, 1),l)])(L 2 )®fe Ck,,j,_, 


we conclude from Proposition 2.5, (13) and (14) the orthogonal decomposition 

00 k 

F = E|F| + E yy E . ...M) ■.. ,3k-,]). 

k=l 1=0 k; jA,_i 


( 15 ) 


□ 
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Remark 2.8. We deduce from (15) that 


c,(F) = E|F] + Y.Y.Y.T. ....• • ■.A-il) 


k = l 1=0 k; 


In order to compute the expectation of products of multiple integrals (see formula (16) 
below) we introduce some notation following [TT], [20], [IZ] and [T5] . 


• If n G N* then [n] := {1, • • • ,n}. 

• For J C [n] we denote by Of the singleton containing that x G {0,1}” for which 
Xi = 0 i E J holds. 


• If ni,..., n; (/ G N*) are given and n := rii + ■ ■ ■ + ri; we will denote by tb the ’natural’ 
partition of [n] given by the summands rij : 


:= {{1,..., Til },..., {ni H-h ni_i + 1,... ,n}}. 


• Let Tin denote the set of all partitions of [n] and 11* denote the set of all subpartitions 
of [n]. 

• Let n(ni,... ,nz) C n„ (respectively n*(ni ,... ,ni) C 11*) denote the set of all a G n„ 
(respectively a G 11*) with Itbi fl J| < 1 for 1 < i < / and all J E a. 

• Let n> 2 (ni ,... ,ni) (respectively n= 2 (ni,..., ni)) denote the set of all a E n(?7,i ,... ,ni) 
with I J| > 2 (respectively \J\ =2) for all J E a. 

• In order to distinguish between integration w.r.t. the Brownian motion and compen¬ 

sated Poisson process we consider for C [n] {J^ will stand for integration w.r.t the 
Brownian motion) and introduce 11=2,ni,... ,ni) as the set of all pairs (r, a) of 
subpartitions from 11* (rii ,... ,ni) such that for all J G r: | J| = 2 and Ujer J = QS 
well as for all J G a: | J| > 2 and Ujeo- = N \ ■ 

• For r G If* let |r| = ff{J C [n] : J G r} i.e. the number of its blocks and ||r|| : = 

#Ujer J- 

• For (r,a) G 'n.= 2 ,> 2 {J^'-,ni,... ,ni) and / : ([0,T] x {0,1})” —)■ M we dehne frua '■ 

[0, —)■ M by identifying the time variables of each block of r U a and setting 

Xj = 0 for z G Ujst j Xj = 1 for z G Ujeo- J- 

Example: Let rzi = 2,rz2 = 2 and rza = 3. Then T = {{1,2}, {3,4}, {5, 6, 7}}. If = 
(2,4, 6, 7} and r = {{2, 6},{4, 7}}, a = (1, 3, 5} we change by rUa the function /((ti, Xi), • • • , 
{tj,X'i)) into 


frVJaiW ^2, h) = f{{ti, 1), (ts, 0), (ti, 1), (tg, 0), (tl, 1), (^2, 0), (ts, 0)). 
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Proposition 2.9. Let G ((5o + K(5i)) (ui G N for 1 < i < 1) be symmetric and 

assume that for all (r, a) G n= 2 ,> 2 (</^; ■ ■ •, it holds 


Then 



< cx). 

rUcr 


EnU/n.(/nJ= E E 


' (8)/n (iAl"l + l‘"l. 


Proof. Let us assume for the moment that the /„. are of the form 


(16) 


- = Ti-kL^di{tk,Xk) 


(17) 


for some di G L^(A ® (5o + ^f^i))- 

For = ffJi we let J^o denote the multiple integral of order w.r.t. the Brownian 

'‘'i 

motion and I^i the multiple integral of order nj {nj := Ui — n^) w.r.t. the compound Poisson 
process. Similar to (12) we get 



^ i 

Consequently, since 

E = xG{0,1}"% 


JrQrii] 


E EnLi/n,((®fclld.)lo„-J 


JlC[ni],...,JiC[n;] 




JlC[ni],...,J;C[n;] 


i: Ein'.,/.5(idi(.,o)]®’>?)]Ein'../r.(K(',i)]®"‘)|. 

' I I 

JlC[ni],...,J;C[n;] 

From mi Corollary 7.3.2] we conclude 







while m Theorem 3.1] (see also |2ni Section 3.2]) implies 





So we have shown relation (16) for the special situation 10 where each /„. is given as 
tensor product df^\ The general assertion follows by approximation using the multilinear 
nature of ^ w.r.t. (/„,,..., /nJ- □ 
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2.2. Hermite and Charlier polynomials 
2.2.1. Hermite polynomials 

Let us introduce the Hermite polynomials {Km)meN defined by 

^ t, X e M. 

m>0 


With the convention iL_i = 0 we have the relations K'^[x) — xK'^{x) + mKm{x) = 0 
and K!^{x) = Km-i{x), for all m G N. The normalized sequence {\/^.Km)m&n forms an 
orthonormal basis in L^(M, p), where /i denotes the normalized centered Gaussian measure. 

Every square integrable random variable F, measurable with respect to , admits the 
following orthogonal decomposition 

= + Kn. i f e.(s)dB,) . (18) 

fc.i |„|.t i>i V-'« / 


where n = (nj)j>i is a sequence of non-negative integers, |?7,| := Y.iyiF and (ej)j>i is 
an orthonormal basis of L^(0,T). Taking into account the normalization of the Hermite 
polynomials we use, we get 


do = E [F], dl = n\¥. 


F X 


IliM 



s)dBs 


where n • — riixi (nj). 


Now we choose G N and let {toHi,--- ,fv} be a regular grid of [0,T], i.e. Vz G 
{0,..., N}, ti = ill where h = ^. From now on we will use a fixed orthonormal basis (ej)j>i 
of T^(0, T) : we set 

ei{t) ■■= zG{0,...,iV} (19) 

and complete this sequence to a basis in L^(0,T), for example, by using the Haar basis on 
each interval Let = (nf,...,n^) be the vector of non-negative integers such 

that \n^\ = k. Then (see m Proposition 5.1.3]) 





o 



n 


B\ N 


\n 


B\ 




I * 2 = 1 



where ABi = Bj. — Bi._^ and o stands for the symmetric tensor product. 


( 20 ) 


2.2.2. Charlier polynomials 

Definition 2.10. The Charlier polynomial of order m G N and of parameter t > 0 is defined 
by 


Go(x, t) = 1, Ci{x,t) = X — t, X G M 


and by the relation 


C'm+i(x,f) = (x - m - t)Cm{x,t) - mtCm-i{x,t). 
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The sequence 


. ^ Cm{-, Kt)) is an orthonormal basis for L^{N, where 

denotes the law of a Poisson random variable with parameter nt. Let = (nf,..., n^) be 
the vector of non-negative integers such that \n^\ = k. Using the same grid and the same 
functions (ej)i<i<Ar as for (20), we have (see m Proposition 6.2.9]) 


N 




(ei 1 o ... oe^ ^) = - 


n^\ \ h' 2 i=i 


HQ- (AiV„«h) 


( 21 ) 


where AAj = Nj. — The following Lemma gives some useful properties of the chaos 

decomposition. 

Lemma 2.11. 

• Let F be a r.v. in Vp > 1, we have E(|Cp(F)p) < E(|Fp). 


• Let H be in Hf.(M). ITe deduce from Remark 2.8 that Cp Hgds^ = Jq Cp{Hs)ds. 

• For all F G for all i G {0,1} and for all t < r, Df^E,[Cp(F)] = E,[Cp_i(Df^F)]. 
2.3. Truncation of the basis 

Instead of summing over all G and jfe_; G we only consider the N hrst functions 
(ei,...,eAr) of the basis (e*)* dehned in (19). This gives (together with the orthogonal 
projection onto the chaos up to order p) the following approximation of F 

C^(F) = E|F] 


k=i 1=0 kie{i,- ,Ny 

( 22 ) 

Let us now rewrite C^{F) {p < N) in terms of Hermite and Charlier polynomials. From 
(111, (@, d^) and ([2l|) we derive using the notation of Lemma |2.7| that 


{9k, e[(/ci, 0),..., (fcz, 0)] ® e[(ji, 1),..., {jk-i, l)])(L 2 )®fc 

,E (fIIK^b {Gi)C^p{Q,,nh)], 


n 


B\ 


|n|!(Kh)l-^l/2 


i=l 


where we used Gi := and Qi := AAj. From Lemma 2.7 we get then 


N 


C^{F) = do + E E n Knf (Gi) C^b(Q,, Kh) 

k=l \n\=k i=l 


(23) 


where do = E(F) and 


d7 := 


n 




N 


nP\{Kh)\^^\ 


E FnF> (Gi)C„p(<3i,Kft) 


(24) 


2=1 
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Proposition 2.12. Let F be a real random variable in \j‘^{Ft) and let r be an integer in 
{1, • • ■ , N}. For all U-i < t <tr, we have 


Ei (C^f) = do+ 

^ ''t — tr-l 


E E d 

k=l |n(r)|=/c 


~Y' 


KnB 


B,-B 


tr—1 




CnP{Nt - Nj^_,,n{t - tr-l)) 


r—1 


X (l[K^B{G,)C^piQ,,nh) 


\i<r 


:=Ar 


b!“>e, (c~(f)) 


i: d 




r—1 


k=l |n(r)|=/c 
n^>0 


h 


K, 


B,-B 


tr — \ 


n^-l 


t — t 


CnP{Nt - Ni^_„K{t - tr-l))Ar 


r—1 


Di‘>E, (C«(F)) 


r—1 


^ ft-t 

E E 

k=l \n{r)\=k 

nf >0 


h 


K„ 


Bt - B,^ 




n^CnP-l{Nt - Ni^_^,K{t - tr-l))Ar 


r—1 


where for r < N n{r) = (n^(r), n^{r)), and n^{r) stands for {nf ,..., nf), where A = B or 
P and Hr = {n^,nf). 

Proof. The first result comes from [3 Proposition 2.7] for the Brownian part and from the 
fact that E,t{Cn{Qr, kih)) = Et[/„(l|’"' ^ ^ j)] = Cn{Nt — K{t — tr-i)) (see [IS Proposition 
6.2.9]). The second result comes from [HI Proposition 2.7]. To get the last one, we write 
Di^^Cr,p{Nt - -tr-i) = (see US Definition 

6.4.1]). 

□ 


Remark 2.13. For t = tr and r >1, Proposition 2.12 leads to 


E,„ (C^F) = do + x: <nA'„f (G,)%(Qi,«:A) 

k=l |n(r)|=fc i<r 

(C"f) = h-'P j2 Y. (Gr) C„f(Qr, Kh) (n A-„f (G.) C„,.(Qi, Kh) 

k=l \n{r)\=k \i<r / 

nf>0 

D!%. (Ct'F) “EE (Go) ntGne-l(Qr, Kfc) f n Anf (G,) C„,(<3„ kA)) 


k=l \n{r)\=k 

nP >0 


\i<r 
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When r = Q, we get (C^ P'j = do and we define (Cp f'^ = (which is 

the limit of {Cp f'^ when t tends to 0) and (Cp F^ = where ei : = 

(1, 0, • • • ,0) of size N and On is the vector null of size N. 


The following Lemma, similar to Lemma 2.11, gives some useful properties of the operator 




Lemma 2.14. Let F be a r.v. in W{F'j’) and H he in 


Then 


. V(p,iV) e (N*)^ E{\C^{F)\^) < E(|C,(F)|2) < E(|Fp), 

• Let H be in Hf{R). We deduce from ( |2^ that Cp (^Jo Hgds^ = Jo Cp {Hs)ds. 

• For all F G for all i G {0,1} and for all t < r, Df^E,[C^(F)] = ¥.r[C^_fiDf F)]. 
Let us end this subsection by some examples. 


Example 2.15 (Case p = 2). From (23)-(24), we have 


N N 

C^{F) =do + E {dJ^KfiGj) + dj^CfiQj, Kh)) + 5] + di^^GfiQ,, nh)) 

j=i i=i 

N j-1 

+ E E {d^^’''KfiGi)KfiGfi + «h)C'i(Q,-, nh)) 

j=li=l 
N N 

+ E E d^^’^’^KfiG,)GiiQj, Kh) 

i=l j=l 


where 


dE = e(fxi(g,)), dJf = ^E(FC'i(g,-,Kh)), 

dE = 2E(FiF2(G,)), dE = 

= ¥.{FKfiGf)KfiGfi), = ^E(FC'i(Q„«:h)C'i(Q„/s:h)), 

= ^nFKi{Gi)GfiQ,,Kh)). 

Remark 12.131 leads to 

^tfiC^{F)) =do + E {dj^FfiGj) + dJ^GfiQj, nh)) + ^ (dE^2(E) + d^^GfiQ,, nh)) 

i=i i=i 

+ EE(4’""i^i {G,)KfiGj) + di^^FcfiAN,, Kh)GfiQj, nh)) 

j=ii=i 

+ E E «:h). 

i=l j=l 


16 









3. Numerical scheme 


3.1. Picard’s approximation 

Picard’s iterations: (F°, U^) = (0, 0, 0) and for q eN, 

e + r / (s, F/, Zl Ut) ds - r ZI+HB, - f Ur^dNs, 0<t<T. 

Jt Jt J]t,T] 

It is well-known that the seqnence Z'^, IP^) converges exponentially fast towards the 
solntion {Y,Z,U) to BSDE Q. We write this Picard scheme in a forward way. Let 
denote := ^ + f (s, F/, F|, f/|) ds. We dehne 


y^^+i = E ( F*? 


a)-J‘ f(s,Y,\ZI,Unds, 


Z?^^ = E f 


Ft- ) , f/r' = E ( 




Ft-]. 


(25) 

(26) 


3.2. Chaos approximation 

Let (F'^’P, denote the approximation of (F”?, Z*?, U^) bnilt at step q using a 

chaos decomposition up to order p\ (F°’P, Z°’^, U^’P) = (0,0, 0) and 


y<?+l,P ^ E 


7^+1,P ^ ^ 


Cp(F'?’^ 


Ft 


Fr^4(F'?’P) 


Ft- 


f f (5, F/’^ zr, Ur) ds, 

Jo 


jj g+hP ^ ^ 


Ftr 


(27) 

(28) 


where F'^’P = ^ + Jq f (s, F/’^, ZfP, Uj'P) ds. 


3.2.1. Truncation of the basis 

The third type of approximation comes from the truncation of the orthonormal L^(0,T) 


basis (ej)j>i dehned in (19). Instead of considering the whole basis we only keep the hrst 

->N 

-'p 

gives us explicit formulas for Et(C^F), F|°^Ef(C^F) and F{^^Ei(C^F). From ([2^ and 


N functions (ei, • • • , cn) to build the chaos decomposition projections . Proposition 


2.12 


28), 


we build {Y^^’PF, Z^'P’^, U^'P’^)g in the following way : (F°’P’^, Z^P’^, f/PpF) = (g, 0, 0) and 


F 


q+l,p,N _ 


E^(C^(F9>pF)) _ f' f Yf’P’^, ZfP'^, ds, 


^q+i,P,N ^ Di^\Et{C^{F‘^’P’^))), = Di^\Et{C^{F'^’P’^))) 


(29) 

(30) 


where := ^ + f(s, Yf’P’^, ZfP’^, Uf’P'^)ds. 

It is not necessary here to use predictable projections of Z^^^'P’^ and fact, 

Z9+1.PF and t/9+PpF are adapted and cadlag, and from their explicit representation given 
above one concludes that the predictable projections are the left-continuous modihcations: 
E^_Zg+i,PF _ ^ 9 _+i,P,iv ^^_^g+i,p,iv _ ^9 _+i,pF_ Moreover, the integral in (|^ does not 

change if one uses left-continuous modifications. 
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3.2.2. Monte Carlo approximation 

Let F denote a r.v. of In practise, when we are not able to compnte exactly 

do and/or the coefficients of the chaos decomposition (23)-(24) of F, we nse Monte- 
Carlo simnlations to approximate them. Let (F"*)i<m<M be a M i.i.d. sample of F and 

be a M i.i.d. sample of ((Gi, Qi), • • ■ , (Gat, Qn))- 
We approximate the expectations of (24) by empirical means 


1 M _ 


n 


B\ 


M 


N 


m=l 




i: F” !!*■«• (CD C„p(Qr.'=A) 


(31) 


2=1 


In the following, we denote 


c"'"(F) = S + y y S” n A'„.(Gi)c„p(o„Kft). 

k=l \n\=k l<i<N 


(32) 


Ei(Cp’“(F)) and Fi(Ei(Cp’“(F))) denote the conditional expectations obtained in Propo¬ 


sition 


2.12 


when (do,c?fc)i<fc<p,|n|=fe) are replaced by {do, dl)i<k<p,\n\=k) ■ 


E, (Cf “f) = Io+ 

E E 

t'l\n^\=k V h 




Kn? 


B,-Bj 


F^^Ei (C;'’^(F) 


= h-^/^Y. E ? 

k=l \n(r)\=k 
ntt>0 




r—1 


h 


K, 


iy tip _ 


'C„p(lV,-iV,,_.,K(i-A_i)) 

X fn''XGi)G„p(0.,K/i)), 


\i<r 


— Ar 


Bt-B 


tr—1 


n^-l 




CnP{Nt - Nj^_,,l^{t - tr-l))Ar, 


r—1 


(C^'^{F)) 


k=l |n(r)|=A: 
nt! >0 


t tj. _"I 

F 




Bt-B 


^r —1 




n^C^P_l{Nt- Nj^_^,K{t-tr-l))Ar. 


r—1 


Remark 3.1. As pointed out in Remark 3.2], when M samples of {F) are needed, 
we can either use the same samples as the ones used to compute do and d'] or use new 
ones. In the first case, we only require M samples of F and (Gi, • • ■ ,Gn,Qi, ■ ■ ■ ,Qn)- 
The coefficients and do are not independent o/ni<i<Ar KnB{Gi)Cnp{Qi, n.h). In this case. 
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the notation E,t{Cp’^{F)) introduced above cannot he linked to "E {C^’^F\iFi^. In the sec¬ 
ond case, the coefficients and do are independent of K^B{Gi)C^p{Qi, nh) and we 

have Et F\Fi^. This second approach requires 2M samples of F and 

(Gi, • • ■ , Gn, Qi, • ■ ■ , Qn)- Convergence results are proved when using the second approach. 

We introduce the processes ^ uq+^,p,N,M'^^ useful in the following. It 

corresponds to the approximation of when we use Cp’^ instead 

of Cp , i.e. when we use a Monte Carlo procedure to compute the coefficients d^. 


yq+l,p,N,M ^ _ j f (gq,P,N,M^ 

yq+l,p,N ,M _ jj{0) ^pq,p,N,M^',^',^ jjq+l,p,N,M _ jj{l) ^pq,p,N,M'^'j-<^ 


where F^-pF.a^ ;= ^ + fj' f(efP’^’^)ds and = (s, ^q,p,N. 


\M 


(33) 

(34) 


the solution of ([^ — and (y 


q,p,N,M 


4. Convergence resnlts 

We aim at bounding the error between (y, Z) 

dehned by (33)-(34). Before stating the main result of the paper, we introduce 
some hypotheses. 

Hypothesis 4.1 (Hypothesis 'Hm)- Let m E N*. We say that F satisfies Hypothesis Hm if 
F satisfies the two following hypotheses 

• : Vj eW F E V^’G i.e. ||F||W < cx) 

• .' Vj E N* V/o,^i e N such that l = lo + li + I<rn there exist two positive 
constants and kf (j) such that for all multi-indices a = (oi, ■ ■ ■ ,azo) ^ 

7 = ( 71 , • • • 5 77 + 1 ) ^ {0, and for a.e. U E [0, T], Si E [0, T] it holds 


ess sup ess sup E|F)" _ 

^i+1 )■■■ + 


F - D'^ 

s*+Zi 


Fffi < kf{j) \ti - s, 


\jhF 


In the following, we denote Kf^{j) = maxz<m kf (j). 

Remark 4.2. If F satisfies Hm, for all I < m and for all multi-indices a = (oi, • • • ,«/) E 
{0,1}^ we have for a.e. (ti, ■ ■ ■ ,ti) E [0,T]^ and (si, ■ ■ ■ ,si) E [0,T]^ that 

|E(£>“....,,F) - E(D“.....„F)| < At(l)(|ii - + ■ ■ ■ + |(, - s/'). (35) 

Hypothesis 4.3 (Hypothesis Hpp). Let {p, N) E We say that a r.v. F satisfies Hpp 

if 


Vpp{F) := V(F) 


EE 

k=l |n|=A: 


(n 


S'] 


N 


m 


{Khffi 


rY{FY[K^BiGGG^piQi,Kh)] <oo, 


2=1 


where Y{ffi) denotes the variance of a r.v. G 
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Remark 4.4. If F is hounded by K, we get Vp^j^{F) < K‘^YX=q (^k \ Hence every bounded 
r.v. satisfies Hp j^. 

This Remark ensues from E (llili HlB{Gi)Clp{Qi, nh)^ = ' • 


Remark 4.5. Let X be the M.-valued solution of 


Xt = x+ / b{s,Xs)ds+ / a{s,Xs)dBs + / 7 ( 5 ,tG[0,T], 

Jo Jo Jo 

where b,a,'j : [0,T] x M —)■ M are (7°’™ functions, uniformly Lipschitz w.r.t. x and Holder 
continuous of parameter | w.r.t. t, with linear growth in x and with bounded derivatives. 
Then every random variable ^ of type piX^) or g Xgds^ with g G C“(M) satisfies Hm 
and Hpj^f for all p and N. (To prove that Urn is satisfied one can use ITR Theorem 3], while 
Hp jv i^ implied by Remark 4-4 ) 

Theorem 4.6. Let m he an integer s.t. 1 <m <p + l. Assume that f satisfies Tip+q+i o.nd 
HIj, and f e C' 0 ,p+g+i,p+g+i,p+g+i^ 


II _ y<1’PX,M y _ jj _ jjg,p,Af,M^||2 


L 2 


^0 

< — + 
- 29 


Ai{q,m) 


(p + 2 - m) ■ ■ ■ (p + 1 ) 


+ 242(g,p) + 


, Afiq,p,N) 


M 


where Aq is given in Section 4H_- Ai is given in Proposition 4F, A 2 is given in Proposition 
4 . 13 , and A^ is given in Proposition \4.15[ 

If f ^ ^0,00,00,00 ^ satisfies "Hoc Cind Hfiro,oo’ SH 

lim lim lim lim ||(F - ^ jj _ jjq,p,N,M\\i 2 ^ g_ 

g^-oo p^.oo AT^-oo M^.oo " ' 


Proof of Theorem \4-(^ We split the error into 4 terms : 

1. Picard’s iterations : 8 ^ = \\{Y — Y'^, Z — U — U^)\W 2 , where (W?, Z^, f/9) is dehned 
by 

2 . the truncation of the chaos decomposition : = ||(^y' 9 _y 9 ,p^ _ 2'9 —t/ 9 _[/ 9 ,p^|| 2 ^^ 

where (Y^’P, Z'I’P, U^’P) is dehned by 

3. the truncation of the 1^(0, T) basis : E'^’PH = ||(F 9 ,p _ ^ z< 1’P - ZfiP'^,,W'P - 

UfP’^)\\l 2 , where (yI’PX, Z^’PH is dehned by 

4. the Monte-Carlo approximation to compute the expectations : 8 '^’PHX = ||(^y<?,p,A^ _ 

yq,p,N,M z‘^’PX _ z'^'PXXi jjq,p,N _ jjq,pXM\\\2 -^pg^e (Y'^’PHXi Z'^'PXX jjq,pXM'\ js 

dehned by (g-Q. 

We have 

II _ y<i,pXM z — zP’PHH jj _ ||2^ ^ 4(£!9 _j_ gQX _|_ gg,pX _j_ gq,pXM^ 


It remains to combine (36), Proposition 4.9, Proposition 4.13 and Proposition 4.15 to get 
the hrst result. □ 
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4.I. Picard’s iterations 

The first type of error has already been studied in 1211 (see the proof of Lemma 2.4), we 
only recall the main result. 

From pn Lemma 2.4], we know that under Hypothesis the sequence (y?, Z"?, U^)q 
dehned by (25)-(26) converges to (Y,Z,U) dP x dt a.e. and in Sf 
Moreover, we have 


X 


X 


||(r -Y^Z- Z<, U - U‘>)\\h < 


(36) 


where Aq depends on T, ||^|p and on ||/(•, 0, 0, 0)11^2 


(0,T) 


4-2. Error due to the truncation of the chaos decomposition 

We assume that the integrals are computed exactly, as well as the expectations. The 
error is only due to the truncation of the chaos decomposition Cp introduced in ([^. 

For the sequel, we also need the following Lemmas. We postpone their proofs to the 
Appendix [Appendix A. 1 


Lemma 4.7. Let m G N*. Assume that f satisfies T-Lln+q ^ £ 0 ,m+q,m+q,m+q ^ Then 

Wq' < q, Wp E N, Z'^' , IT^') and Z‘>'’P, belong to . Moreover, 

IKF", Z'', [/'^)||W < C{j, ||e|L_ . mpf\Uk<m+q). 

Lemma 4.8. Assume that f satisfies TLp and f G Then it holds for any 

J>1 

< C{p,J, m\p,l,{\\d^pf\\oo)k<n.Vp). 

Proposition 4.9. Let 1 < m < p+1. Assume that^ satisfies PLln+q ^ ( 70 -"^+'?-™+'?-”*+'?^ 

We recall 8^'^ = ||(H^ — Z^ — Z^'^, IT^ — We get 


gq+i,p < c'iT(T + l)LlS'^’P + 


A'i(g,m) 


(p + 2 - m) ■ ■ ■ (p + 1) 


(37) 


where Cl is a scalar and Ki(q,m) depends onT, m, ||^|L , ^{m+q-iy. and on (\\d^„f\\oo)i<k<m+a- 


Since 8^'^ = 0, we deduce from (37) that 8'^’^ < where Ai{q,m) : = 

^ X maxi<i<q iFi(/, m). Then, (W’'?, t/'^’?’) converges to (F^, th^) when 

p tends to 00 in ||(•,•,•)||L 2 (see (|^ for the definition of the norm). 

Remark 4.10. ITe deduce from Proposition ^ ■ 9\ that for allT and Lf, we have limp^ac8'^’^ = 
0. When CiTiT + l)Lj < 1, i.e. forT small enough, and iff satisfies PL]^ and f G ^ 

we also get limp^oolimq^oc8'^'^ = 0. Indeed, it holds lim^^oo < i-cIt(!t+i)l) (p+ 2 -J)-(p+i) 
and supj Ki{j, m) < 00 since from the proof of Proposition 4-9 one concludes that Kfij, m) = 

60(||?iri,». + r/„’'||/(s,n’,ziVi)lli,...<is) <c(r,m, ||^L+„tf±t^,(||S5),/|U<„+j)). 

' (m —1)! 
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Proof of Proposition\4-!^ In the following, we denote := Yff'^ — AZ^’^ := — Zf, 

AUr := - U! and Af^ ■= f{t, Y^, . U^) - f{t, Y,\ Zl Uf). Firstly, we deal with 

E[supo<t< 7 ’From (25) and (27) we get 


=Ei[C„(F'?’P) - F^] - [' Aff’Pds, 

Jo 


=Ea 4 ( 0 -e]+Et 

- [' Af^ds. 

Jo 


C 


r /(». rr,zr,ur)dB\ - jf Ks, y;, zj, t/j)* 


We introduce ±Cp /(s, Yf, Z|, Uf)ds"j in the second conditional expectation. This leads 
to 


AYf^^^'P =Ea4(0-e]+Et 


C 


Afrds 


Af^ds 


E/ 


4 (/(5, Yf, Zl Uf)) - f{s, Yf, Zl Uf)ds 


where we have used the second property of Lemma |2.11| to rewrite the third term on the 

r.h.s. 

From the previous equation, we bound E[supo<i<r \AYl^^’'^\l by using Doob’s maximal 
inequality and the Lipschitz property of / 


II sup |Ay,'^+'’"|||2<2||C,(0-el|2 + 2 
o<t<r 


C 


Afl^ds 


+ 2 / \\C,{f{s,Yl,ZlUl) - fis^Yf^ZlUlhds 
Jo 

+ Lf r |||AF/’^>| + |AZ,^’^| + |Af/,^’?’||| 2 ds. 

Jo 

To bound the second term on the r.h.s. of the previous inequality, we use the hrst property 
of Lemma 2.11 and the Lipschitz property of /. Then, we bring together this term with the 
last one to get 


sup 

0<t<T 


|ak,''+‘''’||| 2 < 2||C,«) - ^||2 + 2 F ||c,(/(s,y/, z ;, t/’)) - /(s,r;,z’, 

J 0 
T 

+ 3Lj [ |||AF/’?'| + \AZll + lAUllhds. (38) 

Jo 


Let us now upper bound || AZI+^’P^ds + k, Jq || A7/|+^’^’|||(is. To do so, we use the Ito 
isometry || AZ^+i’P^ds + k \\AUf+^’P\\lds = || AZ^+^'PdB^ + AUf+^’PdNs\\l 

Using the Dehnitions (26)-(28) of Uf^^) and (Zf~^^’^, and the Clark-Ocone 
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Formula (see na Theorem 1.8]) leads to 

r AZ^^+^’PdBs + r = F^- E(F'^) - (Cp(F''’^’) - E(Cp(F''’^’))), 

Jo Jo 

= ir ‘ + r /(s. U’. zi ut)ds - y?+‘ 

Jo 

- 1 k ?+'->’ + jF(s,Yr,Zl’‘,Ut-’’)da-YS*^-J^ . 

Rearranging this summation makes appear We get 


f \\AZ^/^’P\\lds + K f ||Af/,^+^'P||2ds < 4|| sup |Ayj''+^’P|||2 
Jo Jo 0 <t<T 

+ 2L} \\AY,^^2 + + WAU^Vs^ . 

(39) 

Since (fj' {{AY^^’Py + IIA^s^’^lh + WAUP’Pydsf < + 1)SP’P, by computing 5 x 


(38) + (39) we obtain 


f+y" < 60||C,«) - f 11" + 60T / ||4(/(s,y;, Z’, U^) - f(a, YJ, ZJ,(7;)||ps 

'j 0 

+ i37 ^(^ + ^) y(y + 1)l)£P'P. 


K 


Since ^ and /(s, R/, Z|, Uj) belong to (^ satishes "Hm+g, f E C'l 
(y?, ZP, IZ^) G 5™’“ (see Lemma 4.7)), Lemma [2^ gives 


0 ,mH-q',m+g,m+g 


and 


£;<?+i.p < 


6 oiieii 


+ 


60T 


(p + 2 - m) ■ ■ ■ (p + 1) (p + 2 - m) ■ ■ ■ (p + 1) Jo 
411(1 + k), 


\\f{s,Y,^,ZlUml^ds 


K 


-T{T + l)L)£^'P. 


Since ||/(s, Y^, Z^, U'i)\\]jUs is bounded by C{T, m, i\\d^pf\\oo)k<m, || (>"^ f^‘^)||^ 2 m) 

(see (A.l), in the proof of Lemma 4.7), Lemma 4.7 gives the result. 

□ 

4-3. Error due to the truncation of the basis 

Fix N E N* and put h = F. Use {po,Pi} = {l{o}) as orthonormal basis of 

L^({0,1}, ^0 + «^^i) and £x an orthonormal basis {ek)km* for -^^([0; ^])-^([0) ^j);-^) 

such that ti = ih for z = 0,1,..., A and 


e* = 1 <^<A. 
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Lemma 4.11. Assume F = E[F] + In{gn) ^ satisfies (35) with m = p. Then 


E\{C: -C,)iF)\^<Kfi\ 


T 

N 


2I3f 


,T 




2=1 


T 

N 


2I3f 


T(l + T)e^ 


where := (with Kfi := Kf{l) from (g;. 

We refer to Section Appendix A.2 for a proof of Lemma |4.11[ 


Lemma 4.12. Assume f satisfies Tip (i.e. Hypothesis f.l) and f G (7°’^’^’^. Then, for 
all integers q > 0, Iq^p := Jq f{s,Yfi’P, ZfiP,UfiP)ds satisfies Tip so that by Remark 4.2 
for all 1 < r < p and multi-indices i^ G {0,1}'' and for a.e. (ti, • • • fir) ^ [0,T]'' and 
(si, • • • , Sr) £ [0, T]'’ we have 


\E{D\fi..^,Jq^p) - E{Dlfi..^Jq,p)\ < + • • ■ + |fr - 

where = | A /3^, and Kfi’P depends on Kf, ||^||p,i, T and on (ll^6p/lloo)l<A:<p- 
We refer to Appendix A.3 for the proof of Lemma 4.12 


Proposition 4.13. Assume that f satisfies Tip and f G We recall ■= IKP^’P — 

Y<],p,n^ Zi^p - U'i’P - \\l 2 . We get 

^g+l,P,N < ^ l'jl2Z^,P,N ^ K2{q,p) (40) 

where is a scalar and K 2 {q,p) depends on Kp, T, ||^||p,i and on (ll^sp/l|oo)l<fc<p- 
Since £^'PF = we deduce from (40) that £ 1 'PF < A 2 {q,p) where A 2 {q,p) : = 

K 2 (g, p)T{T+l)e'^ ^ ^ t/F’^) converges to {Y^^p, Z^'P,U^^p) 

IK') ■) OIIl^ when N tends to oo. 


Proof of Proposition 4-13 In the following, we denote 

•= — Yfi’^ 


:= Zffi^ - Zfifi := 


and 


:= fit, Yr^\ Zffi\ Uffi’^) - fit, Yfifi Zfifi Ur). 
Firstly, we deal with || supo<t<T \AYfi^^’^’^\\\ 2 . From (|^ and (|^ we get 

^Yfi+^’PF ^ Ei[C^(F''’P’^) - CpiF'^’P)] + f AffiP'^ds. 
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By using th e sec ond property of Lemma 2.14, by following the same steps as in the proof of 
Proposition 4.9 and by introducing ±C^(Jq f{s, U^’^)ds), one gets 


sup < 2|ICK) - CJi)h + 2 


0<t<T 


^ Af^ds 

(C" - CJ vr, zr, U»)ds 

+ Lf r |||AF/’P’^| + \AZfP^^\ + lAU^’P'^lhds. 
Jo 


It remains to apply the first property of Lemma 2.14 to get 


sup < 2||C,^(e) - 4(0II 2 + 2 

0<t<T 




^-C) 

p ^p) 


{f{s,Yr,Zr,Ur)ds 


(41) 


+ 34 r IIIAFO'^’^I + |A4^’P’^| + \AU^'P'^\\\2ds. 

Jo 

Let us now upper bound \\AZ^~^^’P'^\\lds + k \\AU^~^^’P’^\\lds. 

Following the same steps as in the proof of Proposition |4.9| , one gets 

r WAZt^’P’^Wlds + K r WAUr^’P’^Wlds 

Jo Jo 

< 4|| sup lAYt^+^’P’^lWl + 2L] ( r ||AF/’P ’^||2 + WAZ^^’P’^y + ||Af/f^’^Usds ) 
0<t<T ■’ \Jo J 


(42) 


Adding 5 x (41) and (42) gives 


^5+i,P,A^<60||(C^-4)(Oil2 + 60 


K - Cs) (0(/(s. y.", zr, uy)di 


K •’ 


Since ^ and Iq^p satisfy (35) (see Remark 4.2 and Lemma 4.12), Lemma 4.11 


gives 


rjn . 2/3^Al 


g<i+i,P,N < gQ ^ ^ ^ 


A ^2 I ( l^Iq,p\2\ I _|_ 




and (40) follows. 


□ 


4-4- Error due to the Monte-Carlo approximation 

We are now interested in bounding the error between (F^.p.^^ 2 <hP,N^ jj'hp,^'^ dehned by 
(g-(g and ^g,p,7v,M^ defined by (Q-(0. is defined by (Q 

and (32). In this Section, we assume that the coefficients are independent of the vector 


(Gi, • • ■ ^Gn), which corresponds to the second approach proposed in Remark 3.1 


Before giving an upper bound for the error, we recall the following Lemma, which mea¬ 
sures the error between and Cp’^ for a r.v. satisfying Eip j^ (see Hypothesis 4.3). 
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Lemma 4.14. Let F he a r.v. satisfying Hypothesis We have 


Moreover, we have E(|C^’^(F)n < E{\F\‘^) + HVp^n{F). 

We refer to Section [Appendix A.4| for the proof of the Lemma. 


Proposition 4.15. Let^ satisfy Hypothesis Hp j^ and f be a bounded function. Let ■ = 

llj'yg.p.AT _ y<1,pHM ^ Zq,pF _ ^ Jjq,P,N _ ip/g 


gq+i,P,N,M Y C^T{T + 1)L)8‘^’PF,m ^ 


Ks{q,p,N) 

M 


where C 3 is a scalar and K 3 {q,p,N) := C 4 (v^,Ar (0 + ^^ll/ilL X]fc=o (^^)) some C 4 > 0 . 
Since £^’PHM = yj^ deduce from the previous inequality that where 

A 3 {q,p,N) := K^iq, p, N) ^ ^ {YP’<iF,m ^ zP’IFM ^uq,P,N,M^ converges to 

{yi^pH^ Z‘I’PF^ Uq,pF'j 


in 


\i 2 when M tends to 00 . 


The proof of Proposition 4.15 is the same as the proof of |Sl Proposition 4.17], except 


that we consider jumps. The jump part is treated as in (42). 


5. Implementation 

5.1. Pseudo-code of the Algorithm 

In this section, we describe in detail the algorithm. We aim at computing M trajectories 
of an approximation of {Y,Z,Lf) on the grid T = {t* = i'^,i = 0, ■ ■ ■ ,A^}. Starting from 
('yO,p,iv,M^ ^o,p,7V.M^ f/0,p,7V,M^ ^ (0,0,0), (Q-(0 enable to get {yq^pHH ^ z^'P^H ^ijq^pFM-^ 

for each Picard’s iteration q on T. In practice, we discretize the integral /q / 

which leads to approximated values of {Y^’PFH^ z’^’PFH^ jjq^pFM'^ computed on a grid. Let 

us introduce (r^, , ,17^. )i<j< 7 v, dehned by 


^yO,p,NM^ y 0 ,P,N,M^ jjO,P,N,M^ ^ 


and for all g > 0 


^ ii 


= -hY,f {t„ V 


T x^q^pPM -ryq,pF ytI’PHH 


7 

5 


i=i 


.ir-YA 


-^q+l,p,N,M 

jq+l,p,N,M 
ii. 


-iN.M /^q^pHM\ 


ti I H 


where := ^ + h E*=i fiU, 

in the algorithm. 


(43) 

Here are the notations we use 
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• q: index of Picard’s iteration 

• Kit: number of Picard’s iterations 

• M: number of Monte-Carlo samples 

• N: number of time steps used for the discretization of Y and Z 

• p: order of the chaos decomposition 

• G Af7v+i,M(I^) represents M paths of computed on the grid T. 

• Z”? G A4N+i,Mi^) (resp. G AfAr+i,M(lR)) represents M paths of (resp. 

jjq,p,N,M^ computed on the grid T. 

Since ^ G ^ can be written as a measurable function of {Bt, Nt)t<T- Then, 

one gets one sample of ^ from one sample of ((Gi, Qi), ■ ■ ■ , {GnjQn)) (where Gi represents 

— " and Qi represents Nj. — Nj. ). 


Algorithm 1 Iterative algorithm 


1: Pick at random N x M values of standard Gaussian r.v., stored in G, and N x M values 
of Poisson r.v. of parameter nh stored in Q. 

Using G and Q, compute (^™)o<m<Ar-i- 
= 0, ZO = 0, = 0. 

for g = 0 : Kit — 1 do 
for m = 0 : M — 1 do 

Compute (F-?)™ = ^ + hY.ti fih m, (Z«),,,„, (G'?),,™) 

end for 

Compute the vector d = (do, (dk)i<k<p,|n|=k) of the chaos decomposition of 


9 

rF — d_ pg) 

m Jn _ 

7 — 

n^\ 

\ M 



10 

for j = 1 : N do 






11 

for m = 0 : M 

— 1 do 





12 

Compute (IE 

{Df> 

'Et 

. Zj 

(C^’^F-?)))™ 

13 

= 


Ei 

=1 f{ti, 


14 



iCKMp^))) 

m 



15 

= 



m 



16 

end for 






17 

end for 






18 

end for 






19 

Return (Y'^'*)o,: = do, 


1 JeijON 

and (G^« 

)o,: = 


^N,ei 
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5.2. Numerical Examples 
5.2.1. First example 

The following example is borrowed from [T3]. We consider a Poisson process N with 
K = 1 and the following BSDE 

dYt = —cUtdt + ZtdBt + Ut{dNt — dt), 

^ = Nt. 

The explicit solution is given by 

(Fi, Zt, Ut) = (W + (1 + c)(T - f), 0,1). 

Figure represents the evolution of with respect to M when 

q = 5, p = 2 and N = 20. 


Convergence of YO w.r.t. M 



Figure 1: Evolution of {ys.p.N.M^ ^q,p,N,M^ respect to M when p = 2, N = 20, q = 5, c = 0.5, 

T= 1 


Table [^gives the computational time needed by the algorithm with this choice for q, p, N 
and for different values of M. We notice from Figure[^that the value of (Y^’P’^N^ 
is close to the true solution from M = 2 x 10^. When M = 2 x 10^, the CPU time is about 
1 minute, which is quite small. 
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M 

10^ 

5 X 10^ 

10^ 

5 X 10^ 

10^ 

2 X 10^ 

5 X 10^ 

10^^ 

CPU time (in s) 

0.253 

1.277 

2.567 

13.24 

26.81 

56.91 

142.75 

283.65 


Table 1: CPU time w.r.t. M when p = 2, N = 20, q = 5, c = 0.5, T = 1 


5.2.2. Second example 

We consider now the following BSDE 

dYt = —[aYt + (dZt + 'yUt)dt + Z^dB^ + UtdNt, 

^ = exp(aT + bBp + cNp). 

The explicit solntion is given by 

= gar+6Bt+cVfg(a+tt±^^^)(r-t)+(e"-l)(K+7)(T-t) 

Z* = E,-{D^Y^) = bY,-, Ut = EUDlYt) = {e^ - 

We choose a = (3 = 0.3, 7 = 0.2, a = —0.1, b = 0.1, c = 0.2, k = 3 and T = 2. For 
these valnes, we get (Fq, f^o) = (6.599,0.66,1.4612). For M = 4 x 10^, p = 2, iV = 50 
and q = 10, we get = (6.560, 0.56,1.294h We plot one path of 

{Y^’N'M, {Zr^'^, Z,),<T and , f/*),<T in Fignres § | § with M = 4 x 10^ 

p = 2, iV = 50 and q = 10. 



Figure 2: One path of ,Y) 
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one path of (Zapp.Zth) 



Figure 3: One path of , Z) 


one path of (Uapp,Uth) 



Figure 4: One path of ,U) 
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Appendix A. Technical results 


Appendix A.l. Proof of Lemmas 4- 7 and 4-8 


Appendix A. 1.1. Proof of Lemma \4 . 7 

Let ii := max(ti, • • ■ , t;). First, we prove by induction that Vg' < g, U^') belongs 

to i.e. Vj > 2 

\\{Y‘'',Z<',U^’)\\lj= Y, i:esssupjE[ sup \Dll..j,Yf \>] + LEl\Dll_.,,Z;’\i]dr 

l<l<m i; I ii<r<T 

+ CnDl[,..,,U;'\^dr] < ao. 


Let r > ti. Using (25) gives 

where := (s. F/-'. (//-i). 

Jti 

Using the Dehnition of and applying Doob’s inequality leads to 

Ehup^|ci;,^^.,,,y;’'PI < cu) ('E||c;;,..^,,,fp| + e , 

where C{j) is a generic constant depending also on T. Analyzing the outcome of the repeated 
Malliavin derivative where for f {6f~^) the chain rule holds while 


= /(5, Uf'-^ + - /(0f-') 

(see, for example, [HI Lemma 3.2]), one can see that the term \Dl[ ... ^^f{6j'~^)\ is bounded 
by a sum of terms of type 


^0 + ^l+^2 

E Kf\\< 

k=l 






where G {0, l}h are vectors of size Ij and h + h + h Y 1. Then, Holder’s inequality gives 


and 


^ ^ esssupE[ sup 

l<Z<mi;e{0,l}* ti<r<T 


(A.2) 


< 


CU) lieiik. + E E ii^Vllio ■ 


1=1 \k=i 
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From (|2§, we get + Sj" Then 






dr 


Using (A.l) yields 


^ ^ esssup [ E[\Dl[...^t^Zf\^]dr 

l<Z<miie{0,l}* Jti 


m //+1 


< 


c(j) (KII4+1.J + E (E KfWL 1 z*'-'. 


1=1 \k=i 


The same type of result holds for {^dr. Combining this result with (A.2) 

gives 




m+1 




\A:=1 


1=1 


Iterating this inequality yields the result. 


Appendix A. 1.2. Proof of Lemma 4■ 8 


We prove it by induction on q. Let r >ti := max(ti, • • • , t;) and := (s, Uf^). 


From (27) we get that 

J ti 


Jti 


where we have used Lemma 2.11 to get the second equality. Applying Doob’s maximal 
inequality leads to 


E[ sup \Dl...,^Yr^^^\^] <C(j) 

ii<r<T V 

+ e( r\Di...,j{erWds]], 


(A.3) 


where C{j) is a generic constant depending also on T. Let us hrst deal with the hrst 


term of the r.h.s. of (A.3), we assume I < p. Following Proposition 2.5, we know that 


p^q,p ^ x;n=o^n(5'n). Then 

OO 

Dl[,.. ^ n(n - 1) ■ ■ ■ (n - / + zi, • • • , zi)), 


n=l 


32 












with Zk = (4,4) and 


^ n(n - 1) ■ ■ ■ (n - / + l)4_z(^„(*, Zi, • • • , zi)), 


n=l 


P-'- (n + /V 

= ^1 ' 4(^n+/(*, ^1, • • • , ^0)- 


71=0 


n! 


Let us denote ^ni(*) ■= S'ni+K*, ''' , From Proposition 2.9 we get 

(ni + Z)! {nj + l)\ 


p-i 

= ^ E (4^ (gn^ ) ■ ■ ■ Inj {guj ) ) 

Til ,••• ,nj=0 

_ ^ (ni + /)! (nj + /)!.^ 


(A.4) 


rii'. 


Tljl 


ni,---,nj=0 ^ 1 - 


nn 


E 




kl 


J®e[n] (T,cr)en^ 2 ,> 2 (<r®;ni,...,nj) 


/ dA^I+H 


< i: EE'-'ErFriiis-i'ii- E 


711,••• ,nj=0 


rii’ 


i=l 


^kl ykl+kl 


Thanks to the assumptions on / and and induction hypothesis, we have F'^’P G D^’^. 
Then, gives that gni+i{zi,--- .z^+i) = lkni+i||oo < 

ni<p-l,we get Eni+zlU < Then 


p-i 


E|i 4 _,(r>l‘,,...,,,p’nit < EE E kHth+i-i 

,7ij=o jsg[n] (T,(T)gn=2,>2(".^®;>^ivi%) 


<4(p,j)(||F’'^’1Ip.i)T 

We have ||F^’P||p,i = Ez<pEi,e{o,i}'esssup^^^...E(|T))y..where 

E(|E(,...,,F'^’^|) < e(|E(,...,EI) +E( r \Dl...,J{er)\ds). 


(A.5) 


By using (|g, we get E (/? | A‘‘.,.-,i,/(»f’’)l*) < C (eUi I| 3 ,VIIco) IKV''’'. Z"’”, U-mii- 
Then 


l|F'''’’lle,i < ll^llel + 1:0 y || 5 ‘y||» 

l<p \A:=1 / 

< c(p,3) I ii5iij,i+x:c (x: ii8j,/iij„) uy-”,z^’’‘,u'>--)n‘i, 

\ 1<P \k=l } 


(A.6) 

(A.7) 


Let us now deal with the second term of the r.h.s. of (A.3). By using (A.l), we get 

E(/V;;,....<,/(e’')i'4 <c'li:iisj,/iii„) ii(v’^’'.z7-p,p’niiL- (a.8) 
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Combining ( |A.5[ ), ( |A.7| ), ( |A.8[ ) and ( |A.3[ ) yields 
E[ sup 

ii<r<T 


<c(p,j) (||^||Ju + Ee EII3.VlPco ll(i'’’^E’’^c/nllgd loss) 

l<p \k=l / 


+ e i: Iisj fii' 


spj 11 CO 




\k=l 


and 


y esssupE[ sup 

l<Z<mi;g{0,l}* ii<r<T 

/ mVp / ^ \ \ 

< C(pj) [ll?IPp,i + E c IEII A/ll~ ) 11(5"’'". Z’’". f^’'’’)llb) ■ (A.9) 


From (28) we get 




Then 


E||A‘‘.,.-,,,Zr''’’PI<ir <c{l E[|C,_,_,(A‘;:.':;;;;;"’F'''’')PI<ir 1 


Using (A.5) and (A.7) leads to 


^ ^ ess sup [ E[\Dl[... 


I 


< C{p,j) ||?||j,i + EC E II A/lli. 11(5'’’’'. z-'”,ui'nw'Y ■ 


1=1 \k=i 
rT 1 


The same type of result holds for ... Combining these results with 

(|A.9l) gives 


/ mVp / I \ 

ll( 5 «+‘.’’, Z’+‘.>’, C/’+‘'’’)|P™., < C(pJ) + E c EII A/llt 11 ( 5 ’’’'. Z’’’. v-nw'Y 




Iterating this inequality yields the result. 


i=i \k=i 


Appendix A.2. Proof of Lemma 4-11 


We will prove the assertion by induction in p G N. Since {C^){F) = (Co) (A) Lemma 4.11 
holds for p = 0. Assume that for p G N* 


p—1 / T' \ p—1 rpi 
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Since 


(C^ - C,){F) = {Cl, - C,_i)(F) + {Pi - P,)(F), 
it suffices to show that 

9 2Pf T^p 

We have Pp{F) = Ip{gp) where we will assume that Qp is symmetric. It holds 

= 4(0 

with 


9p ((^Ipl)) ■ ■ ■ ) (^pPp)) ^ ^ (fi'p(('Pl)5 ■ ■ ■ ) ("Pp))! ^[^1; ■■•5 ^p])l2([0,T]^) 

kpe{i,-,vp 

X e[/ci, kp]{ti, ■ ■ ■ ,tp). 

Then is constant w.r.t. {ti, ■ ■ ■ ,tp) & Ak^ := x ■ • • x A^^ with Aj :=]tj_i, A] since 
e[ki,...,kp] = h“2lAkp- We have by ([^, (10) and (35) that 


E|(F;^-P,)(F)|^ 

= nipigl)-ipigp)\" 

P'-\\h 2 (c/p((-Ai), •••,(•, ip)), e[ki,kp\)Lmo,T]p) - gp{{-, h), v))llL(AkJ 

kp ip 

= EE-'S! \k ^ 5'p(('Si, A), ■ ■ ■ ,{^p,ip)) gp{{', ii),''' , {', ip))dsi...dsp 

V i "'^k„ 

ivp Ip y 

< ^ ( I \gp{{si,ii), - ■ ■ ,{sp,ip)) - gp{{t,,ii), - ■ ■ ,{tp,ip))\dsi...dsp]dti - ■ -dtp 


A2(Akp) 
\ 2 


Akp 


< / ([ ^^'(1^1 --h |tp - Sp|^^)dsi...dsp^ dt,---dtp 




Appendix A. 3. Proof of Lemma f.l2 
We will show that if 


{Y,’^, ^Uf’^) satisfies Pp for a.e. t E [0,T] (A.10) 

(with ~ ^ /^s) then also Iq^p = fl f{s,Y^’P, Z^’P,UfP)ds does satisfy Pp. As /o,p 

is constant, it satisfies Pp. For g > 1 we will use the notation 

- Dlffifl^F) and prove that for 1 < r < p 
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(that "Hp holds for Iq^p is clear). Setting ts_j := max(ti, • • • • • • , s^) and : = 

17r)wehave 


1,-Vs,-Vts_,- 


+ 


/ tiVSjVts— 

rti\/Si\/tS-i 

I {tiASi)\/ts_i 




From the proof of Lemma 4.7 we know that \D^ 
snm of terms of type 


a{i.i jg bonnded by a 


y \Kf\u I icsn’-i icf.‘zri 

^i=i 

where kj G {0,1}^^ are vectors of size Ij with lo + h + h ^ r and tj are snb vectors of 
{fi,..., fj_i, ti A Si, Sj+i,..., Sr}. Holder’s ineqnality and Lemma 4.8 give 


E 


rt-iW Si\/ts_i 


' {tiASi)VtS-i 


D 


o(l:i—1) 




< C(pj, llfllp.l.(l|9i/l|oo)l!<p)|*. - SiL 


For the hrst term on the r.h.s. of (A.11) we notice that 

f 

JtiWsiWts—i 




is bonnded by a snm of terms of type 


/ t S_j \ • -I 

1 \j = l 


H IKp/IU Ift'S-rl \D^:ArDi‘T’/\du 


where 4/^’^, T^} = Z^’P, U^’P}, and kj G {0, l}b are vectors of size /q, h, h + h + 

1 < r while the tj denote the appropriate snb vectors of {fi,..., fi_i, U, Si, Sj+i,..., s^}. 

By Holder’s ineqnality and assnmption (A. 10) we conclnde that 


E 


' tiVsiVts_i 




< Kr{j)\ti - Sip^A.P. 


We hnish the proof of Lemma 4.12 by arguing that assumption (A. 10) holding for true 
for a certain q, implies it for g + 1 : We want to use (27) and (28) and therefore we hrst 
notice that in the same way as above for Iq^p one can show that (A. 10) implies that 


f f (s, Zl'P, UI'P) ds satishes V,p. 

0 

It is also clear that satisfying Up is stable with respect to linear combination and taking the 
conditional expectation E^. What we still need to check is whether satisfying Tip is also stable 
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with respect to the truncation Cp. For this, let us assume that F = In{9n) satisfies 'Hp. 


Following the proof of Lemma |4.8| , we have 

^OLi P 

CO 

= ^ u(u - 1 ) ■ ■ ■ (n - r + l)In-r{,gn{.*, Zi, ■ ■ ■ , Zi, ^'+ 1 , • ■ ■ , 4 ) 

n=r 

- 9n{*,Zi,--- , 4 )), 


where Zj = {tj,ij) and z' = Like in (A.4) we get 

<C(p,,.T) 'f + 


ni,-" ,n,=0 


ni'. 


Tljl 


X 


n \\9ni+r{*, ^1, • • • , Zi, 4+1, • • • , 4) - ^ni+r(*, ^1, ' ' ' , 4 ' ' ' ) 4 


r// II CO 


2 = 1 


<c(p,i,r)(A'ni)|(j-Sj|'*o' 


where we used that 


(ni + r)!|4„,+4*,Zi,-- - ,Zi,4+i,-- - ,4)-^„.+4*,zi,-- - ,Zi_i,4--- ,4 


r// II CO 


< A-ni)|ii-s. 


I^F 


Appendix A.4- Proof of Lemma 4-14 


Using the definitions (23) and (32) leads to 

p 


N 


_ + _ + + ^ ^ K„pGi)C„pQ„Kh). 

k=l 1721=/;; 


2=1 


Since is independent of (Gj, (5i)i<i<v 

E(|(C" - C"’")(F)|=) = E(|4 - *1^) + 'tY, 

k=l \n\=k 




{n^)\ 


-E(K-d4|4 


The definition of the coefficients do given in (24) leads to 


E(|(C~ - c“)(F)r) = vpo) + E E 

fc=l|n|=fc )■ 

Using the definition of d^ (see (31)) leads to the first result. To get the second result, we 
write C^^^{F) = - C^){F) + C^{F). Since E - C^){F)C^{F)) = 0, we get 

E(|Cp^’^(F)n = E(|(Cp^’^ -Cp^)(F)n +E(|Cp^(F)n. 


Lemma 2.14 ends the proof. 
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Appendix A.5. The product of two multiple integrals 


For the convenience of the reader, we cite here [T^ Theorem 3.6] from Lee & Shih adapted 
to our simple situation where the multiple integrals Ik{gk) are built using the process B + N 
like in Subsection 2.1.2 For this, we hrst introduce the ’contraction and identification 
operator’ For symmetric functions pk G (L^)®^(A (g) (ho + ^hi)) and Pm G 0 

(ho + «^hi)) we define the function pk 
([0,T]x{0,l})'^Mby 


thru : ([0,T] X {0,1})“ X ([0,T] x {0,1}) 


m—a—b 


X 


( 5 'fc ®afi'm)(x,y,z) = /" 5 (fc(x,z,w) 5 (m(w,z,y)d[A 0 (ho + Khi)]®“(w) (A.12) 

2([0,T]x{0,1})“ 

for (x,y,z) G ([0,T] x {0,1})^““-^ x ([0,T] x {0,1})™-““^ x ([0,T] x {0,1})^ 

Theorem Appendix A.l. If Pk ^ (L^)®*’(A 0 (ho +/thi)) and pm ^ (L^)®’”(A 0 (ho + fthi)) 
are symmetric functions such that \pk\ 0^ \pm\ is in (^/, 2 ^( 8 i(fc+m- 2 a-b)j'y ^ then 


kAmkAm—a /L^ 

Ikigk)Imigm) = 

a=0 b=0 \^/ 





^k-\-m—2a—b {gk ®a g^ 


An immediate consequence is that if pk and pm have disjoint support, then Ik{pk)Im{gm) = 
Ik+m{gk 0 gm) ■ 
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